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Abstract 

We present a new family of gauge invariant non-local order parameters 
for (non-abelian) discrete gauge theories on a Euclidean lattice, which are 
in one-to-one correspondence with the excitation spectrum that follows from 
the representation theory of the quantum double D{H) of the finite group 
H. These combine magnetic flux-sector labeled by a conjugacy class with an 
electric representation of the centralizer subgroup that commutes with the 
flux. In particular cases like the trivial class for magnetic flux, or the trivial 
irrep for electric charge, these order parameters reduce to the familiar Wilson 
and the 't Hooft operators respectively. It is pointed out that these novel 
operators are crucial for probing the phase structure of a class of discrete 
lattice models we define, using Monte Carlo simulations. 
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1. Discrete gauge theories and their excitations 

In two-dimensional gauge theories we distinguish two classes of particle- 
like excitations, which we call electric and magnetic. Electric excitations 
are labeled by nontrivial representations of the gauge group, and are either 
put into the theory as external charges or as dynamical fields in the action. 
The magnetic excitations are labeled by topological quantum numbers, on 
the classical level related to solitonic sectors of the gauge theory. In planar 
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physics these are the magnetic fluxes and non-abehan generahzations thereof. 
On the quantum level these manifest themselves as particle-like excitations 
carrying these topological quantum numbers. 

The magnetic excitations may also carry electric charges, which leads to 
exotic particles called dyons. In the two dimensional setting these particle 
types may behave like anyons with fractional spin obeying highly nontrivial 
braid statistics. Generally speaking, quantum groups provide the language in 
which the quantum physics of two-dimensional anyonic systems is optimally 
casted and understood. In that sense a complete classiflcation of all sectors 
of a discrete gauge theory (DGT) can be achieved within the mathematical 
framework of the quantum double D{H) [1] of the discrete gauge group H 

Experimental realizations of discrete gauge theories are still modest, though 
there are interesting proposals for implementing them in Josephson junction 
networks 0, 0, [sl and in systems of polar molecules in optical lattices [gj. 
Globally symmetric implementations were explored in [?, [8|. In a larger 
context it should be pointed out that DGT's play a crucial role in Kitaev's 
[sl seminal paper on topological quantum computation, as they constitue 
the most elementary examples of Topological Quantum Field Theories. This 
simplicity follows from the fact that there is no underlying Conformal Field 
Theory involved as is the case for Chern-Simons theories. 

1.1. Transformations on DGT states: quantum symmetry 

We will not give a detailed account on the emergence of quantum group 
symmetry in DGT, this can be found in the literature but do present 



a short summary of the basics to flx the notation and introduce some key 
concepts required later on. 

Consider the following operators acting on states in the Hilbert space of 
a DGT. First there is the flux projection operator, denoted by Ph, which acts 
as follows on a state I ■?/'): 



if the state contains flux h 
otherwise 



Secondly, we have the operator g, for each group element g E H, which 
realizes a global gauge transformation by the element g: 

9 m = \y) , 



2 



where it should be noted that we have not yet modded out by the gauge 
group to obtain the physical Hilbert space. 

These operators do not commute, and realize the following algebra: 

PhPh' = Sh,h'Ph 
gPh = Pghg-^g ■ 

The set of combined flux projections and gauge transformations {Phg}h,geH 
generates the quantum double D{H), which is a particular type of algebra 
called a Hopf algebra. In the following we will see that its representations 
correspond one to one with the electric, magnetic and dyonic sectors of the 
discrete gauge theory. 

1.2. Representation theory of the quantum double: particle spectrum 

The representation theory of the quantum double D{H) of a finite group 



H was first worked out in but here we follow the discussion presented in 
lo| and follow the conventions of that paper. 

Let ^ be a conjugacy class in H. We will label the elements within A as 

{%,%,■■■ G A, (2) 

for a class A of order k. In general, the centralizers for the different group 
elements within a conjugacy class are different, but they are isomorphic to 
one another. Let "^A^ G H he the centralizer for the first group element in 
the conjugacy class A, ^hi. 

The set relates the different group elements within a conjugacy class 
to the first: 

= l^^/j^j^T/jj, ■ ■ ■ j'^/ij. I "^hi = "^Chi "^hi ^h!^} ■ (3) 

This still leaves a lot of freedom, but we fix our convention such that "^Xh^ = e, 
the group identity element. The centralizer ^N, being a group, will have 
different irreps, which we label by a. The vector space for a representation 
a is spanned by a basis "vj. The total Hilbert space that combines magnetic 
and electric degrees of freedom, V^, is then spanned by the set of vectors 

iWrv,)} , (4) 

where i runs over the elements of the conjugacy class, i = 1,2, ■ ■ ■ , dim A and 
j runs over the basis vectors of the carrier space of a, j = 1,2, ■ ■ ■ , dim a. 
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To see that this basis is a natural one to act on with our flux measurements 
and gauge transformations, consider an irreducible representation 11;^ of some 
combined projection and gauge transformation Phg: 

m 

where the element g is the part of the gauge transformation g that commutes 
with the flux ^/ii, defined as 

9=\lg-^9'^Xh,- (6) 
This indeed commutes with the element ^hi. 

''hi~g=^h^\lg-.g^xh^ 

=^Xgh,g-i g % g'^ g ^xh^ 



Note that when performing a successive series of gauge transformations 
the definition ([6]) of the centralizer element working on the charge part en- 
sures that the left 4c element of the first transformation equals the right 4c 
element of the second, so that one is left with terms like 

(^a;(gg/)ft^(g<;/)-l g' ^ gh,^ ^ ^^^^ gh,g 9 h,) ^. 

(7) 



(^^i99')hr{.gg')-^ ^9' 9) ^Xh, 



mj 



which is a property required later on in the definition of the lattice order 
parameter. 

2. The lattice formulation of a DGT 

We will first briefly introduce the spacetime lattice approach to gauge the- 
ories. This approach to gauge field theories was introduced in the context of 
quark confinement [12|], but some research concerning discrete gauge groups 



131 . |14| has been conducted too. We note however that the latter papers 
were unable to capture the full richness of the quantum double spectrum. In 
particular, the dyonic sector of the theory was not touched upon. 
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For a more exhaustive introduction to lattice gauge theory, consider the 



review articles [15|, [ly, [ITj or textbooks [18|, [19 . 

After setting the stage we introduce a new class of order parameters 
corresponding to any of the sectors in the quantum double representation 
theory. 

2.1. Lattices, gauge fields and actions 

We discretize spacetime into a set of sites i, j, ■ ■ ■ using a rectangular 
lattice. The gauge field Uij, which takes values in the gauge group H, lives 
on the links ij, jk, ■ ■ ■ connecting sets of neighboring sites. The links are 
oriented in the sense that Uij = . 

We note that the gauge field Uij takes care of the parallel transport from 
site i to site j, were we to have included a matter field charged under the 
gauge group. An ordered product of links along a closed loop measures the 
holonomy of the gauge connection. 

Gauge transformations are labeled by a group element gi E H and are 
performed at the sites of the lattice. The gauge field transforms as 

Uij ^ gi Uij gj^ . (8) 

The standard form for the lattice gauge field action makes use of the or- 
dered product of links around the smallest closed loop possible, the plaquette 
ijkl: 

Up = Uijki = Uij Ujk Uki Uii , (9) 
which transforms as a conjugation under gauge transformations 

Up^g,Upgr\ (10) 



The action per plaquette, which corrsponds to the well-known form F 



in the continuum limit for H = SUn [12], is given by 



S=-J2Pe.xAUp) , (11) 

a 

where Xa is the group character in irrep a and Pa is inversely proportional 
to the coupling constant for irrep a. For SUn gauge theories one usually 
only includes the fundamental representation and is thus left with only one 
coupling constant. However, we will not make this restriction here. Gauge 
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(a) Construction of the operator A^. (b) Location of the operator A^: set of pla- 

quettes S corresponding to a closed loop on 
the dual lattice, with the associated loop on 
the real lattice. Also shown is the site io 
acting as the basepoint. 

Figure 1: Two building blocks of the operator A^ and the full operator. 



invariance of the action is ensured by the fact the characters are class func- 
tions, in other words we consider actions where the number of independent 
couphngs equals the number of classes i.e. the number of irreps. 

With these definitions, the (imaginary time) path integral becomes 

Z = I DUe-^v ^(^p) = n E n ^"''^'''^ • (12) 

2.2. The order parameters 

In gauge theories, the order parameters are non-local operators in one- 
to-one correspondence with the different particle species. Gauge invariance 
of these operators is a fundamental requirement, and furthermore Elitzur's 
theorem [20| guarantees that the quantum expectation value of a non-gauge 
invariant operator vanishes. 

Here we present a class of order parameters labeled by a class A and 
centralizer irrep a that have a one-to-one correspondence with the particle 
species in the spectrum of the quantum double representation theory. 
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2.2.1. Construction 

We will first show what the order parameter locally looks like and after 
that give a definition of the operator as a whole. It is well-known that in 
the lattice picture a magnetic flux corresponds to a nontrivial group element 
living on a plaquette. On the other hand, an electric charge consists of a 
lattice site that lies on a closed loop of links, of which the traced product 



is evaluated in the path integral [l9j. In this paper we construct a loop 



operator combining both flux and charge, taking proper care of the intricate 
interlinkage between them ([5]) we alluded to before. 



Consider Figure 1(a), We have marked the site Zq, which will act as a 



reference point or flux bureau of standards [21| . The sites Iq to 12 are where 
the charge lives. Two links are named Uqi and U12 and the two plaquette 
products are Up^ to Up^. We take the plaquette product to start at the point 
where the vertical links touch the plaquette and multiply the edges of the 
plaquette clockwise. Note we have left out the plaquette at the bottom; our 
building block for the operator will be the link pointing upward followed by 
the plaquette. 

We now want to insert fluxes labeled by a class A into the plaquettes pi 
and p2 and make sure that a charge a will live at the sites io, h and 12- To 
insert a flux, we have to 'twist' the path integral. What we mean by this is 
that we want to make sure that in the integration over field configurations 
the dominant contribution will not come from the configuration where all 
the plaquettes take the trivial value in the group, but from the one where all 
plaquettes are trivial except the chosen plaquettes pi and p2 (note that we 
will need to extend this set of plaquettes to one corresponding to a closed 
loop on the lattice; the present discussion serves as an intermediate step to 
arrive at the final definition.) 

Let us pick a flux state h in the conjugacy class A. We will compare 
fluxes at the point ^o- This means that the holonomy around a plaquette pi 
with basepoint zq will equal h for all plaquettes pi. The flux to be inserted at 
plaquette pi is thus U^i hUoi, whereas the flux to be inserted at p2 equals 
U{,'U^,'hUoiU^2. 



We will now start at io in Figure 1(a) and work our way upwards. We 



associate the link Uqi with the plaquette pi. This association is necessary 
because we require a flux state to map an element in the full group H to the 
centralizer "^N. The representation matrix to be put on the first link is 

D{x^^UoiXu-i^uJ , (13) 
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where we have used equation with Uqi playing the role of g, to map 
elements of the full group to the centralizer. 

Next is the insertion of flux at pi. We want the path integral to favor 
a configuration with a nontrivial group element at this plaquette instead of 
the trivial one. The action favors the trivial element, so by changing the Up-^ 
appearing in the action to {U^i hUoi) ^ Up^, the dominant contribution to 
the path integral will come from configurations where Up-^ = Uq^ hUoi. The 
relevant operator expression to put in the path integral is thus 

where the first term eats the original Boltzmann factor appearing in the path 
integral and is replaced by the second. 

Up to the next representation matrix. This one is associated with the 
plaquette p2, where we insert flux Uqi h Uqi Uu- This gives a matrix 

D(X^}_1,^^ U12 XTT-lrr-lhTT rj ). (15) 

Note that the left element between brackets of this expression equals the 
inverse of the right element between brackets in equation f|T3|) . This is crucial 
since it will allow us, when we make use of the representation property of 
these matrices, to keep gauge invariance. 

The flux insertion at plaquette p2 is completely analogous to the former 
case. This leaves us with with the following expression for (a small part of) 
the order parameter: 

■ ■ ■ ^'(-^^;.w.,f^-%r.V,rw.,.J«^'""'-^'''" ■ (16) 

This makes clear the idea of the construction. Now we will introduce some 
extra notation to make the general expression not too convoluted. The dif- 
ferent fluxes in a conjugacy class A are called h^. The correct group element 
to be inserted into the action at plaquette pj is denoted by h^^ . For example, 
in equation (fT4|) . this group element would be U^i Uqi Up-^. 

If we now draw a closed loop on the dual lattice, this loop pierces a 
set of plaquettes. We call this set S. With this notation and the above 
considerations, the anyonic operator is given by: 
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The sum over fluxes within a class is required for gauge invariance, as we will 
show next. 

2.2.2. Gauge invariance 

Gauge invariance of expression fll7l) is assured by two of its properties. 
First, every building block contains a product of group elements taking care 
of the parallel transport to a single reference point io- This already ensures 
gauge invariance with respect to all local gauge transformations, except at 
the point zq. 

Gauge invariance at zq is guaranteed by summing over the group elements 
within the class A. A gauge transformation at this point by an element g 
changes the inserted flux from hi to ghig~^, which corresponds to a reordering 
in the sum over elements in A. 

One might argue that singling out a reference point zq seems artificial. 
However, the same argument that we just gave can be used to show that 
the point we pick is arbitrary. Selecting a neighboring point i'^ would give a 
conjugation of hi by the element Uigi'^, which would result in a reordering of 
the terms in the sum over elements in A. 

2.2.3. Closing the loop 

We built the order parameter by considering a building block of a single 
link and a plaquette. The full order parameter is a closed loop of these 
building blocks, in the sense that the links form a closed loop on the lattice 
and the plaquettes are pierced by a closed loops of links on the dual lattice. 

To go around a corner, it might be necessary to include two links or no 
links at all in a single building block. This is straightforward. However, how 
to close the loop on itself might be less trivial at first sight. 

Since each Da in equation (fT7|l is a representation matrix of the central- 
izer group, we can use the property Da{g)Da{h) = Da{gh). This results 
in all of the Xh^ factors canceling in the definition of the operator, except 
for the first and the last. However, since total group element looks like 

U\oopXjj-i it is still of the form (jS]) and is thus an element of the 

centralizer group. 
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2.2.4- Topological spin 

The dyonic sectors in a DGT posess a feature called topological spin [lO 
This means that under a 27r rotation, their wavefunction obtains a nontrivial 
spin factor (for ordinary bosons or fermions this factor would be +1 or — 1 
respectively, but in a topological theory such as a DGT more general complex 
phases are allowed.) For a dyon carrying flux A and charge a, the topological 
spin is defined as 



XaiA) 



(18) 



da 

where da is the dimension of the irrep a. 

The topological spin appears when the electric part of the loop on which 
the operator is evaluated winds around a single plaquette once. Let us call 
the loop with such a winding S'. The same loop without the winding is called 
H. The topological spin appears as 

(^f(2')) = e^--<--)(Df(S)), (19) 

since the electric charge a picks up a flux in class A in this operator expec- 
tation value. 

2.2.5. Comparison with known order parameters 

The order parameters that so far have been used in lattice gauge the- 



ories are the Wilson loop [12| and, to lesser extent, the 't Hooft loop [22 
The former corresponds to the creation and later annihilation of a particle- 
antiparticle pair, whereas the latter does the same for a fiux-antifiux pair. 
The order parameters proposed here are not simple products of these two op- 
erators, but one can 'disentangle' the flux and charge part by making special 
choices for either the conjugacy class A or the irrep a. 

When we select for the class the trivial class in equation (ITTl) , the expres- 
sion reduces to the expression for the Wilson loop Wa- 

Dl = Wa = {U,,i,Ui,,, ■ ■ ■ , (20) 

where it should be noted that in the case of the trivial class (or any class 
defined by a central group element) the centralizer is the whole group, and 
the set of allowable centralizer representations is given by the full set of 
representations of H. 
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On the other hand, when a trivial group representation is selected in 
equation (ITTI) . the expression reduces to 

D^ = H'' = ^ \2 n e^^^^^^-^^'^'''''^-^ , (21) 
which is an 't Hooft loop for class A. Our definition here is slightly different 



from the one appearing in the literature [13j, where the sum over the class is 
performed independently for each plaquette. The authors in that reference 
perform these sums independently to symmetrize over the different possible 
places in the plaquette where the flux insertion can take place, since e.g. 

h U12U23U34. 1/4,1 ^ Ui2hU23U34.Uii etc. 

Now, since we provide a framing, the flux insertions are all performed at 
the point where the electric loop touches the different plaquettes. One could 
achieve the same effect with deflnition f[T7|) . by performing a number of dif- 
ferent twists at each plaquette (loop Q around the plaquettes in S) and then 
picking the trivial centralizer irrep. However, since only dyons have a non- 
trivial topological spin fio'], this action is immaterial. 

2.2.6. Computational complexity 

The usefulness of formulating this order parameter in the context of Eu- 
clidean lattice gauge theory (instead of e.g. using the Hamiltonian approach 
jol) lies in the fact one can use Monte Carlo (MC) methods to study its 
behavior and probe the phase structure in the parameter space of the theory. 

As equation f|T71) does look quite complex one might fear that the opera- 
tors we have proposed are not very convenient to use in MC simulations. We 
like to argue out that this is not the case and will further substantiate this by 
explicit calculations in a future publication. We flrst like to point out that in 
MC simulations most CPU time is by far used for generating new conflgura- 
tions. The time it takes to perform a measurement of an operator like (fT7I) 
within a single conflguration is negligible compared to that of generating a 
new fleld conflguration. The second point is concerned with obtaining a good 
sampling resolution on the generated MC conflgurations. The problem here 
is that the dominant contribution to the expectation value of the operator 
we propose will not coincide with the dominating conflgurations generated 
by the MC algorithm in the absence of the operator. This problem has been 
studied in the context of f/i gauge theory and magnetic monopole operators 



23|, and their results are directly applicable here. 
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3. Conclusion and outlook 



We have proposed a family of order parameters for discrete gauge 
theories on a Euclidean lattice. They are in one-to-one correspondence with 
the excitation spectrum that follows from the representation theory of the 
quantum double D{H) of the finite group H. Since they explicitly contain 
the magnetic flux label A and electric charge label a, they directly allow 
quantum double physics to be studied on the lattice. In particular cases 
these order parameters reduce to the familiar Wilson and the 't Hooft loop 
operators respectively. 

These order parameters can be studied using Monte Carlo simulations. 
In a forthcoming paper, we will report on the phase diagram of a non-abelian 
discrete gauge theory using the new order parameters proposed here. 

The authors would like to thank Jan Smit and Joost Slingerland for illu- 
minating discussions. JCR is supported by the Stichting voor Fundamenteel 
Onderzoek der Materie (FOM) of the Netherlands. 
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